Mapping

Package Dependencies/ Data

library(sf)
library(ggplot2) #development version!
## devtools::install_github("tidyverse/ggplot2")
library(tidyverse)
library(readr)
library(cowplot)
library(sp)
library(gridExtra)
library(dplyr)
library(ggrepel)
library(plyr)
library(leaflet)
library(rgdal)
library(leaflet)
library(ggmap)

#arrancoordinates <- read.csv("../alldata/ukpostcodes.csv") %>%
# filter(substr(postcode,1,4)=="KA27")

DZBoundaries2016 <- read_sf("../alldata/SG_SIMD_2016")

arran2016 <- DZBoundaries2016[c(4672,4666,4669,4671,4667,4668,4670),]
#Reorder arran 2016
reorderedvector<- c("S01011174", "S01011171", "S01011177", "S01011176", "S01011175", "S01011173", "S01011172" )
arran2016 <- arran2016 %>%
  slice(match(reorderedvector, DataZone))

simple.sf <- st_as_sf(arrancoordinates, coords=c('longitude','latitude'))
st_crs(simple.sf) <- 4326
exampleshapes <- sf:::as_Spatial(arran2016$geometry)
examplepoints <- sf:::as_Spatial(simple.sf$geom)
examplepoints <- spTransform(examplepoints, CRS("+proj=longlat +datum=WGS84"))
exampleshapes <- spTransform(exampleshapes, CRS("+proj=longlat +datum=WGS84"))
namingdzpostcode <- over(exampleshapes, examplepoints, returnList = TRUE)

listID <- list(1,2,3,4,5,6,7)
function100 <- function(argument) 
{
  argument <- arrancoordinates[namingdzpostcode[[argument]],] %>% mutate(DataZone=argument)
}
newarrancoordinates <- lapply(1:7,function100)
newarrancoordinates <- rbind(newarrancoordinates[[1]], newarrancoordinates[[2]], newarrancoordinates[[3]], newarrancoordinates[[4]], newarrancoordinates[[5]], newarrancoordinates[[6]], newarrancoordinates[[7]])
newarrancoordinates$listID <- revalue(as.character(newarrancoordinates$DataZone),
               c('1'="S01004409/S01011174", '2'="S01004372/S01011171", '3'="S01004353/S01011177", '4'="S01004352/S01011176", '5'="S01004351/S01011175", '6'="S01004350/S01011173", '7'="S01004349/S01011172"))

///

Coordinates

postcodelist <- paste(unique(newarrancoordinates$listID), "Postcodes", sep=" ")
datazonelist <- paste(unique(newarrancoordinates$listID), "Datazones", sep=" ")
m = leaflet() %>% addTiles()

m %>% setView(-5.227680, 55.582338, zoom = 10) %>% 
#allcoordinates
addMarkers(
    lng = newarrancoordinates$longitude, lat = newarrancoordinates$latitude,
    label = newarrancoordinates$postcode,
    labelOptions = labelOptions(noHide = F), group = "Postcode Plots") %>%
hideGroup("All Postcode Plots") %>% 
#alldatazones  
addPolygons(data=exampleshapes, 
            weight = 2, 
            label = datazonelist,
            group = "All Datazones") %>% 
hideGroup("Datazones") %>% 
  
#selectcoordinates
addMarkers(
    lng = newarrancoordinates$longitude, lat = newarrancoordinates$latitude,
    label = newarrancoordinates$postcode,
    labelOptions = labelOptions(noHide = F), group = newarrancoordinates$listID) %>% 
hideGroup(newarrancoordinates$listID) %>% 
#selectdatazone
addPolygons(data = exampleshapes[1] , 
            weight = 2, label = datazonelist[1], group = datazonelist[1]) %>% 
addPolygons(data = exampleshapes[2] , 
            weight = 2, label = datazonelist[2], group = datazonelist[2]) %>% 
addPolygons(data = exampleshapes[3] , 
            weight = 2, label = datazonelist[3], group = datazonelist[3]) %>% 
addPolygons(data = exampleshapes[4] , 
            weight = 2, label = datazonelist[4], group = datazonelist[4]) %>% 
addPolygons(data = exampleshapes[5] , 
            weight = 2, label = datazonelist[5], group = datazonelist[5]) %>% 
addPolygons(data = exampleshapes[6] , 
            weight = 2, label = datazonelist[6], group = datazonelist[6]) %>% 
addPolygons(data = exampleshapes[7] , 
            weight = 2, label = datazonelist[7], group = datazonelist[7]) %>% 
hideGroup(datazonelist[1]) %>%
hideGroup(datazonelist[2]) %>%
hideGroup(datazonelist[3]) %>%
hideGroup(datazonelist[4]) %>%
hideGroup(datazonelist[5]) %>%
hideGroup(datazonelist[6]) %>%
hideGroup(datazonelist[7]) %>%
#Layers control
addLayersControl(
    baseGroups = c("All Datazones", "Postcode Plots", "Nothing"),
    overlayGroups = c(newarrancoordinates$listID, datazonelist),
    options = layersControlOptions(collapsed = TRUE)
  )

Inputing example markers.

cliniccoordinates <- read.csv("../alldata/clinics.csv") %>%
dplyr::left_join(arrancoordinates, by="postcode")

Column `postcode` joining factors with different levels, coercing to character vector

#change to character
cliniccoordinates$X <- as.character(cliniccoordinates$X)

m %>% setView(-5.227680, 55.582338, zoom = 10) %>% 
#allcoordinates
addMarkers(
    lng = newarrancoordinates$longitude, lat = newarrancoordinates$latitude,
    label = newarrancoordinates$postcode,
    labelOptions = labelOptions(noHide = F), group = "All Postcode Plots") %>%
hideGroup("All Postcode Plots") %>% 
#alldatazones  
addPolygons(data=exampleshapes, 
            weight = 2, 
            label = datazonelist,
            group = "All Datazones",
            highlightOptions = highlightOptions(color = "black", weight = 2,
      bringToFront = TRUE)) %>% 
hideGroup("All Datazones") %>% 
  
#selectcoordinates
addMarkers(
    lng = newarrancoordinates$longitude, lat = newarrancoordinates$latitude,
    label = newarrancoordinates$postcode,
    labelOptions = labelOptions(noHide = F), group = newarrancoordinates$listID) %>% 
hideGroup(newarrancoordinates$listID) %>% 
#selectdatazone
addPolygons(data = exampleshapes[1] , 
            weight = 2, label = datazonelist[1], group = datazonelist[1]) %>% 
addPolygons(data = exampleshapes[2] , 
            weight = 2, label = datazonelist[2], group = datazonelist[2]) %>% 
addPolygons(data = exampleshapes[3] , 
            weight = 2, label = datazonelist[3], group = datazonelist[3]) %>% 
addPolygons(data = exampleshapes[4] , 
            weight = 2, label = datazonelist[4], group = datazonelist[4]) %>% 
addPolygons(data = exampleshapes[5] , 
            weight = 2, label = datazonelist[5], group = datazonelist[5]) %>% 
addPolygons(data = exampleshapes[6] , 
            weight = 2, label = datazonelist[6], group = datazonelist[6]) %>% 
addPolygons(data = exampleshapes[7] , 
            weight = 2, label = datazonelist[7], group = datazonelist[7]) %>% 
hideGroup(datazonelist[1]) %>%
hideGroup(datazonelist[2]) %>%
hideGroup(datazonelist[3]) %>%
hideGroup(datazonelist[4]) %>%
hideGroup(datazonelist[5]) %>%
hideGroup(datazonelist[6]) %>%
hideGroup(datazonelist[7]) %>%
  
#cliniccoordinates
addMarkers(
    lng = cliniccoordinates$longitude, lat = cliniccoordinates$latitude,
    label = cliniccoordinates$X,
    labelOptions = labelOptions(noHide = F), group = "All GP clinics") %>%
  hideGroup("All GP clinics") %>%   
#cliniccoordinates
addMarkers(
    lng = cliniccoordinates$longitude, lat = cliniccoordinates$latitude,
    label = cliniccoordinates$X,
    labelOptions = labelOptions(noHide = F), group = cliniccoordinates$X) %>%
  hideGroup(cliniccoordinates$X) %>%   
  
#Layers control
addLayersControl(
    baseGroups = c("All Datazones", "All Postcode Plots", "All GP clinics", "Nothing"),
    overlayGroups = c(newarrancoordinates$listID, datazonelist, cliniccoordinates$X),
    options = layersControlOptions(collapsed = TRUE)
  )

pal <- colorNumeric(
  palette = "Blues",
  domain = exampleshapes2$Percentile)
leaflet() %>% 
  addTiles() %>% 
  setView(-5.227680, 55.582338, zoom = 10) %>% 
#alldatazones  
addPolygons(data=exampleshapes, 
            weight = 2, 
            label = datazonelist,
            group = "2004",
            highlightOptions = highlightOptions(color = "black", weight = 2,
      bringToFront = TRUE)) %>% 
hideGroup("2004") %>% 
  
addPolygons(data=exampleshapes, 
            weight = 2, 
            label = datazonelist,
            group = "2006",
            highlightOptions = highlightOptions(color = "black", weight = 2,
      bringToFront = TRUE)) %>% 
hideGroup("2006") %>% 
  
addPolygons(data=exampleshapes, 
            weight = 2, 
            label = datazonelist,
            group = "2009",
            highlightOptions = highlightOptions(color = "black", weight = 2,
      bringToFront = TRUE)) %>% 
hideGroup("2009") %>% 
  
addPolygons(data=exampleshapes, 
            weight = 2, 
            label = datazonelist,
            group = "2012",
            highlightOptions = highlightOptions(color = "black", weight = 2,
      bringToFront = TRUE)) %>% 
hideGroup("2012") %>% 
addPolygons(data=exampleshapes2, 
            weight = 2, 
            label = datazonelist,
            group = "2016",
            highlightOptions = highlightOptions(color = "black", weight = 2,
      bringToFront = TRUE)) %>% 
hideGroup("2016") %>% 
  
#cliniccoordinates
addMarkers(
    lng = cliniccoordinates$longitude, lat = cliniccoordinates$latitude,
    label = cliniccoordinates$X,
    labelOptions = labelOptions(noHide = F), group = cliniccoordinates$X) %>%
  hideGroup(cliniccoordinates$X) %>%   
  
#Layers control
addLayersControl(
    baseGroups = c("2004", "2006", "2009", "2012", "2016", "Nothing"),
    overlayGroups = c(cliniccoordinates$X),
    options = layersControlOptions(collapsed = TRUE)
  )

exampleshapes2 <- as(arran2016, "Spatial") %>%
spTransform(CRS("+proj=longlat +datum=WGS84"))
pal <- colorNumeric(
  palette = "Blues",
  domain = exampleshapes2$Percentile)
leaflet() %>% 
  addTiles() %>% 
  setView(-5.227680, 55.582338, zoom = 10) %>% 
addPolygons(data = exampleshapes2, 
            weight = 2, 
            label = datazonelist,
            group = "2016",
            fillOpacity =0.8,
            color = ~pal(exampleshapes2$Percentile),
            highlightOptions = highlightOptions(color = "black", weight = 2,
      bringToFront = TRUE))

exampleshapes3 <- as(arransimd, "Spatial") %>%
spTransform(CRS("+proj=longlat +datum=WGS84"))
pal2 <- colorNumeric(
  palette = "Blues",
  domain = exampleshapes3$Percentile)
#exampleshapes3$year=2016
#exampleshapes3[exampleshapes3$year == 2016, ]
  
leaflet() %>% 
  addTiles() %>% 
  setView(-5.227680, 55.582338, zoom = 10) %>% 
addPolygons(data = exampleshapes3[exampleshapes3$year == 2016, ], 
            weight = 2, 
            label = datazonelist,
            group = "2016",
            fillOpacity =0.8,
            color = ~pal2(exampleshapes3$Percentile),
            highlightOptions = highlightOptions(color = "black", weight = 2,
      bringToFront = TRUE))

arransimd$listID <- revalue(arransimd$DataZone,
               c("S01004409"="S01004409/S01011174", "S01004372"="S01004372/S01011171", "S01004353"="S01004353/S01011177", "S01004352"="S01004352/S01011176", "S01004351"="S01004351/S01011175", "S01004350"="S01004350/S01011173", "S01004349"="S01004349/S01011172", "S01011174"="S01004409/S01011174", "S01011171"="S01004372/S01011171", "S01011177"="S01004353/S01011177", "S01011176"="S01004352/S01011176", "S01011175"="S01004351/S01011175", "S01011173"="S01004350/S01011173", "S01011172"="S01004349/S01011172"))
#Check you have updated arransimd
exampleshapes3 <- as(arransimd, "Spatial") %>%
spTransform(CRS("+proj=longlat +datum=WGS84"))
pal2 <- colorNumeric(
  palette = "viridis",
  domain = exampleshapes3$Percentile)
#exampleshapes3
listlistlist <- paste(datazonelist, exampleshapes3$Percentile, sep=" ") %>%
paste("%", sep="")
leaflet() %>% 
  addTiles() %>% 
  setView(-5.227680, 55.582338, zoom = 10) %>% 
#alldatazones  
addPolygons(data=exampleshapes3[exampleshapes3$year == 2004, ], 
            weight = 2, 
            label = listlistlist[29:35],
            group = "2004",
            fillOpacity =0.8,
            color = ~pal2(Percentile),
            highlightOptions = highlightOptions(color = "black", weight = 2,
      bringToFront = TRUE)) %>% 
hideGroup("2004") %>% 
  
addPolygons(data=exampleshapes3[exampleshapes3$year == 2006, ], 
            weight = 2, 
            label = listlistlist[22:28],
            group = "2006",
            fillOpacity =0.8,
            color = ~pal2(Percentile),
            highlightOptions = highlightOptions(color = "black", weight = 2,
      bringToFront = TRUE)) %>% 
hideGroup("2006") %>% 
  
addPolygons(data=exampleshapes3[exampleshapes3$year == 2009, ], 
            weight = 2, 
            label = listlistlist[15:21],
            group = "2009",
            fillOpacity =0.8,
            color = ~pal2(Percentile),
            highlightOptions = highlightOptions(color = "black", weight = 2,
      bringToFront = TRUE)) %>% 
hideGroup("2009") %>% 
  
addPolygons(data=exampleshapes3[exampleshapes3$year == 2012, ], 
            weight = 2, 
            label = listlistlist[8:14],
            group = "2012",
            fillOpacity =0.8,
            color = ~pal2(Percentile),
            highlightOptions = highlightOptions(color = "black", weight = 2,
      bringToFront = TRUE)) %>% 
hideGroup("2012") %>% 
addPolygons(data=exampleshapes3[exampleshapes3$year == 2016, ], 
            weight = 2, 
            label = listlistlist[1:7],
            group = "2016",
            fillOpacity =0.8,
            color = ~pal2(Percentile),
            highlightOptions = highlightOptions(color = "black", weight = 2,
      bringToFront = TRUE)) %>% 
hideGroup("2016") %>% 
  
#cliniccoordinates
addMarkers(
    lng = cliniccoordinates$longitude, lat = cliniccoordinates$latitude,
    label = cliniccoordinates$X,
    labelOptions = labelOptions(noHide = F), group = cliniccoordinates$X) %>%
  hideGroup(cliniccoordinates$X) %>%   
addLegend("bottomleft", pal = pal2, values = exampleshapes3$Percentile,
    title = "SIMD Percentile",
    labFormat = labelFormat(suffix = "%"),
    opacity = 1
  )  %>%  
#Layers control
addLayersControl(
    baseGroups = c("2004", "2006", "2009", "2012", "2016", "Nothing"),
    overlayGroups = c(cliniccoordinates$X),
    options = layersControlOptions(collapsed = TRUE)
  )

listlistlist

 [1] "S01004409/S01011174 Datazones 71%" "S01004372/S01011171 Datazones 58%"
 [3] "S01004353/S01011177 Datazones 50%" "S01004352/S01011176 Datazones 79%"
 [5] "S01004351/S01011175 Datazones 56%" "S01004350/S01011173 Datazones 78%"
 [7] "S01004349/S01011172 Datazones 52%" "S01004409/S01011174 Datazones 64%"
 [9] "S01004372/S01011171 Datazones 66%" "S01004353/S01011177 Datazones 78%"
[11] "S01004352/S01011176 Datazones 76%" "S01004351/S01011175 Datazones 52%"
[13] "S01004350/S01011173 Datazones 67%" "S01004349/S01011172 Datazones 50%"
[15] "S01004409/S01011174 Datazones 60%" "S01004372/S01011171 Datazones 63%"
[17] "S01004353/S01011177 Datazones 68%" "S01004352/S01011176 Datazones 63%"
[19] "S01004351/S01011175 Datazones 51%" "S01004350/S01011173 Datazones 61%"
[21] "S01004349/S01011172 Datazones 48%" "S01004409/S01011174 Datazones 59%"
[23] "S01004372/S01011171 Datazones 61%" "S01004353/S01011177 Datazones 71%"
[25] "S01004352/S01011176 Datazones 65%" "S01004351/S01011175 Datazones 51%"
[27] "S01004350/S01011173 Datazones 61%" "S01004349/S01011172 Datazones 54%"
[29] "S01004409/S01011174 Datazones 44%" "S01004372/S01011171 Datazones 56%"
[31] "S01004353/S01011177 Datazones 65%" "S01004352/S01011176 Datazones 63%"
[33] "S01004351/S01011175 Datazones 47%" "S01004350/S01011173 Datazones 61%"
[35] "S01004349/S01011172 Datazones 53%"

Using GoogleMap images instead.

E.g

map <- get_map(location = “texas”, zoom = 6, source = “stamen”) ggmap(map, fullpage = TRUE)

map <- get_map( location = c(-77.0425, 38.8925), # painfully picked by hand source = “google”, zoom = 14, maptype = “satellite” ) ggmap(map)

Can plot over google maps, but not make it interactive.

qmap(location = “arran”)

// //

Trying to make a chloropleth.

bins <- c(0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100) pal <- colorBin(“YlOrRd”, domain = arransimd$Percentile, bins = bins)

arransimd10 <-arransimd st_crs(arransimd10) <- “+proj=longlat +datum=WGS84”

arransimd10 %>% st_set_crs(“+proj=longlat +datum=WGS84”)

addpolygonsarran2016 <- arransimd10 %>% filter(year == 2016, DataZone == “S01011177”)

m %>% setView(-5.227680, 55.582338, zoom = 10) %>% addPolygons(data = addpolygonsarran2016$geometry, weight = 2) %>%

m %>% setView(-5.227680, 55.582338, zoom = 10) %>% addPolygons(data = arransimd10, weight = 2) %>%

addLayersControl( overlayGroups = c(datazonelist, datazonelist[3]), options = layersControlOptions(collapsed = TRUE) )

leaflet() %>% addTiles() %>% addPolygons(filter(arransimd10$geom, DataZone == “S01011177”)) setView(-5.227680, 55.582338, zoom = 10)

addPolygons(data = exampleshapes[1] , weight = 2, label = datazonelist[1], group = datazonelist[1]) %>% addPolygons(data = exampleshapes[2] , weight = 2, label = datazonelist[2], group = datazonelist[2]) %>% addPolygons(data = exampleshapes[3] , weight = 2, label = datazonelist[3], group = datazonelist[3]) %>% addPolygons(data = exampleshapes[4] , weight = 2, label = datazonelist[4], group = datazonelist[4]) %>% addPolygons(data = exampleshapes[5] , weight = 2, label = datazonelist[5], group = datazonelist[5]) %>% addPolygons(data = exampleshapes[6] , weight = 2, label = datazonelist[6], group = datazonelist[6]) %>% addPolygons(data = exampleshapes[7] , weight = 2, label = datazonelist[7], group = datazonelist[7]) %>% hideGroup(datazonelist[1]) %>% hideGroup(datazonelist[2]) %>% hideGroup(datazonelist[3]) %>% hideGroup(datazonelist[4]) %>% hideGroup(datazonelist[5]) %>% hideGroup(datazonelist[6]) %>% hideGroup(datazonelist[7]) %>%

install.packages(“geojsonio”) library(“geojsonio”)

arranDZleaflet <-geojsonio::geojson_read(“../alldata/SG_SIMD_2016/SG_SIMD_2016.shp”, what = “sp”)

leaflet() %>% setView(-5.227680, 55.582338, zoom = 10) %>% addTiles() %>% addPolygons(data = arranDZleaflet, weight = 2)

addPolygons(data = arranDZleaflet, weight = 2, label = arranDZleaflet, group = arranDZleaflet) %>% hideGroup(arranDZleaflet) %>% addLayersControl( baseGroups = c(arranDZleaflet), overlayGroups = c(arranDZleaflet), options = layersControlOptions(collapsed = TRUE) )